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ABSTRACT 

A two-dimensional, elastic-plastic finite-element analysis was used with a critical crack-tip- 
opening angle (CTOA) fracture criterion to model stable crack growth in thin-sheet 2024-T3 
aluminum alloy under monotonic loading after precracking at different cyclic stress levels. 
Tests were conducted on three types of specimens: middle-crack, three-hole-crack and blunt- 
notch tensile specimens. An experimental technique was developed to measure CTOA during 
crack growth initiation and stable tearing using a high-resolution video camera and recorder. 
Crack front shapes were also measured during initiation and stable tearing using a fatigue 
marker-load technique. Three-dimensional elastic-plastic finite-element analyses of these crack 
shapes for stationary cracks were conducted to study the crack-front opening displacements. 

Predicted load against crack extension on middle-crack tension specimens agreed well with test 
results even for large-scale plastic deformations. The analyses were able to predict the effects 
of specimen size and precracking stress history on stable tearing. Predicted load against load- 
line displacements agreed well with test results up to maximum load but the analyses tended to 
overpredict displacements as cracks grew beyond the maximum load under displacement-con- 
trolled conditions. During the initiation phase, the measured CTOA values were high but de- 
creased and remained nearly constant after a small amount of stable tearing. The constant 
value of CTOA agreed well with the calculated value from the finite-element analysis. The 
larger CTOA values measured at the sheet surface during the initiation phase may be associ- 
ated with the crack tunneling observed in the tests. Three-dimensional analyses for nonstraight 
crack fronts predicted much higher displacements near the free surface than in the interior. 

INTRODUCTION 

The aging aircraft research being conducted worldwide is aimed at developing and implement- 
ing advanced fatigue and fracture mechanics concepts into the damage tolerance analysis meth- 


odology for the current and next generation fleets. The phenomenon of "flapping" provided 
the early motivation for much of this research. Flapping is the process by which a small re- 
gion of a cracked fuselage structure peels open, leading to safe decompression. Whereas a 
single crack in a fuselage structure produces flapping, tests on simulated fuselages with multi- 
site or multi-element damage (MSD or MED) are showing fuselage failures without safe de- 
compression [1]. One of the objectives in the NASA Aging Aircraft Research Program [2] is 
to develop the methodology to predict flapping or failure in damaged fuselage structures using 
a finite-element shell code with global-local adaptive mesh capabilities and appropriate local 
fracture criteria to predict progressive failure in complex structures. The present research fo- 
cuses on the development of the local failure criteria using the finite-element method. 

Stable crack growth in metallic materials under mode I (tensile) loading has been studied ex- 
tensively using elastic-plastic finite-element methods [3-10]. These studies were conducted to 
develop efficient techniques to simulate crack extension and to examine various local and 
global fracture criteria. Proposed fracture criteria include crack-tip stress or strain, crack-tip- 
opening displacement or angle, crack-tip force, energy release rates, J-integral, and the tearing 
modulus. Of these, the crack-tip-opening angle (CTOA) or displacement (CTOD) at a speci- 
fied distance from the crack tip was shown to be the most suited for modeling stable crack 
growth and instability during the fracture process [5,9,10]. Newman [11] used critical CTOD 
values obtained from compact specimens to predict failure loads within about 10 percent of 
experimental failure loads for several crack configurations in two aluminum alloys and a very 
ductile steel. Brocks and Yuan [12], Newman et al. [13] and Demofonti and Rizzi [14] found 
that, for various materials and thicknesses, CTOD or CTOA was nearly constant after a small 
amount of crack growth. In some cases, the region of stable tearing where CTOA or CTOD 
was not constant appears to be related to thickness but this region has not been defined quanti- 
tatively and needs to be studied further. 

Numerous investigators [15-20] have also experimentally measured CTOD or CTOA during 
die fracture process. Luxmoore et al. [15] showed that CTOA was constant from the onset of 
stable crack growth in two aluminum alloys, but found different values for different crack 
configurations (middle-crack and double-edge crack tension specimens). Paleebut [16] meas- 
ured CTOD at the initiation of stable tearing in compact specimens made of two aluminum al- 
loys; these results agreed well with numerical values [11] used to model initiation, stable tear- 
ing and instability. Experiments conducted by Schwalbe and Hellmann [17] correlated a 
modified CTOD parameter with crack extension data for various specimen types. Reuter et al. 
[18], using microtopography [19], measured CTOD at the initial crack front location and 
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found a nearly linear relation with crack extension for low-strength steel. These results imply 
that CTO A was nearly constant from initiation. Recently, Kobayashi et al. [20], using a frac- 
ture surface topography analysis (FRASTA), measured CTOA on thin-sheet aluminum alloy 
and found that CTOA was nearly constant after the blunting process. 

The objective of the present paper is to use the critical CTOA fracture criterion to study crack 
initiation and stable tearing in thin-sheet 2024-T3 aluminum alloy material used in aircraft fu- 
selage construction. The fracture sequence in thin-sheet aluminum alloy is illustrated sche- 
matically in Figure 1. At low fatigue stress conditions, the crack front tunnels slightly in the 
interior but the fatigue surface remains relatively flat (indicated by the unshaded area). Then, 
during monotonic loading to failure, the crack front will tear stably under nearly flat fracture 
conditions and begin to develop shear lips (at 45°) through the thickness (indicated by the 
shaded area in Fig. 1). At the end of the transition region, the fracture surface is completely a 
shear mode fracture, tearing on a 45° plane. Under high fatigue stresses, the crack surface 
may have already completed the transition to the shear mode during the fatigue cycling; thus, 
during monotonic loading, the stable crack growth is completely shear mode fracture. 

A two-dimensional (2D), elastic-plastic finite-element analysis (FEA) was used to model the 
fracture process after growing a crack under a precracking stress history. Comparisons are 
made between measured and calculated applied stress against crack extension and applied stress 
against load line displacements for small and large width middle-crack tension M(T) speci- 
mens. A videographic technique was developed to measure the CTOA during initiation and 
stable tearing. The measurements were compared with the critical CTOA values (y c ) deter- 
mined from the 2D finite-element analyses. The effects of precracking stress history on the 
fracture process were also studied. Crack-front shapes were measured during initiation and 
stable tearing using a fatigue marker-load technique. Three-dimensional, elastic-plastic finite- 
element analyses of these crack shapes for stationary cracks were conducted to assess the im- 
pact of non-straight crack fronts on crack-opening displacements. 

MATERIAL, SPECIMENS AND EXPERIMENTAL PROCEDURES 

Tests were conducted on the three specimen types shown in Figure 2, middle-crack, blunt 
notch and three-hole-crack tension specimens. All specimens were made of 2024-T3 alumi- 
num alloy sheet material ranging in thicknesses from 0.05 to 0.09 inches. The 0.09 inch thick 
material was obtained from a special NASA Langley stock (purchased in the 1950's) whereas 
the 0.05 and 0.07 inch thick materials are of a more recent vintage. 
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Fracture tests were conducted on middle-crack tension M(T) specimens (Fig. 2(a)) in all thick- 
nesses. Two specimen half-widths were tested (w = 1.5 and 5.9 inches). The M(T) speci- 
mens were fatigue precracked at either a low or high stress level (at R = 0) to produce an in- 
itial crack. The initial crack-length-to- width (cj/w) ratio was nominally 1/3 for all M(T) 
specimens. The specimens were loaded under displacement control and crack extension (Ac) 
was measured using a 60X microscope during periodic hold times. Crack-opening displace- 
ments were measured at the centerline of the crack. On some M(T) specimens, the crack-front 
shapes were also measured during initiation and stable tearing using a fatigue marker-load 
technique. These tests were stopped after various amounts of crack extension and fatigue (high 
R-ratio) marker loadings were applied. The fatigue precracking and marker loading regions 
resulted in lighter regions on the fracture surface, distinguishable from the darker regions 
which resulted from the stable tearing. The interface between the light and dark regions repre- 
sented the crack front shape after precracking and after stable tearing. 

The blunt-notch specimen, Figure 2(b), was used to verify whether the 2D FEA could model 
large-scale plasticity deformations. These test specimens were similar to the M(T) specimens 
except that a small hole (diameter of 0.05w) was drilled at both ends of a saw-cut (c/w = 
1/3). The holes were polished to help prevent premature fracture. Displacements were 
measured at the centerline and at both notch roots using ring gages. The specimens were 
loaded under displacement control until the specimen fractured. 

A specially designed three-hole-crack tension specimen [21] (Fig. 2(c)) was used to measure 
CTOA in a structurally-configured specimen. The three-hole-crack specimen has a stress-in- 
tensity factor solution like that of a cracked, stiffened panel [22]. The specimen dimensions 
and stress-intensity factor solution are given in Reference 21. The specimen width, w, was 10 
inches and the initial fatigue crack length, Cj, was 0.5 inches. The specimen was loaded under 
displacement control and crack extension (Ac) was, again, measured with the 60X microscope 
during periodic hold times. 

A photographic technique was developed to measure CTOA during crack initiation and stable 
tearing. A high-resolution video camera (500X) connected to a video recorder was used to re- 
cord the fracture process for several middle-crack tension specimens. During crack extension, 
the camera was adjusted to focus ahead of the crack tip so that when the crack grew, it would 
grow into the view of the camera. A typical video frame is shown in Figure 3. In each 
frame, the critical value of CTOA (vj /„) was measured at several distances behind the crack tip. 
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This was done to determine what effect, if any, that varying this distance would have on the 
CTOA measurement. 

FINITE-ELEMENT ANALYSES 

Two elastic-plastic finite-element codes, ZIP2D [23] and ZIP3D [24], were used in the current 
study. The elastic-plastic analysis of both codes employs the initial-stress concept [25] which 
is based on incremental flow theory and small strain assumption. A multi-linear representation 
of the uniaxial stress-strain curve for 2024-T3 was used in the analyses with the von Mises 
yield criterion. The three specimen configurations shown in Figure 2, middle-crack tension, 
blunt notch and three-hole-crack tension, were analyzed. 

ZIP2D, a two-dimensional analysis, was used to study stable crack growth in the M(T) speci- 
mens and large-scale plastic deformations in the blunt-notch specimen under plane-stress and 
plane-strain conditions. ZIP2D uses two-dimensional, constant-strain, triangular elements. 
The mesh pattern used for the 5.9-inch-wide M(T) specimen is shown in Figure 4(a). A simi- 
lar mesh pattern was used for the 1.5-inch-wide specimen. The minimum element size (d) 
along the line of crack extension was the same for all 2D meshes, d = 0.01875 inches. Ficti- 
tious springs were used to change boundary conditions associated with crack extension. For 
free nodes along the crack line, the spring stiffnesses were set equal to zero; for fixed nodes, 
the stiffnesses were assigned extremely large values. See Reference 23 for details of the elas- 
tic-plastic finite element analysis with crack extension. To model fatigue precracking, the in- 
itial crack length was set equal to the initial saw-cut length. Cyclic loads were applied to the 
model and the crack was extended one element length each time the maximum fatigue stress 
(Smax) was reached. Once the crack had grown to the desired crack length, monotonic load- 
ing (displacement control) was applied and the crack growth was governed by a critical CTOA 
criterion. 

A critical CTOA (q> c ) criterion is equivalent to a critical CTOD (5 C ) value at a specified dis- 
tance d behind the crack tip since vj/ c = 2 tan"* (8 c /2d). Whenever the CTOD equaled or 
exceeded a preset critical value (5 C ) during incremental loading, the crack tip node was re- 
leased and the crack advanced to the next node. This process was repeated until crack growth 
became unstable under load control or the crack reached the desired length under displacement 
control. As will be explained later, the critical vj/ c value was analytically determined to match 
the average maximum load measured in several tests of M(T) specimens. 
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The elastic-plastic analyses of M(T) specimens with stationary straight or curved crack fronts 
were done using ZIP3D [24], ZIP3D uses eight-noded, hexahedral elements. Six layers of 
elements were used through the half-thickness (z-direction), as shown in the Figure 4(b). The 
thicknesses of each layer were 0.005615, 0.005635 , 0.01125, 0.01125, 0.00765, and 0.0036 
inches, respectively, where the thinnest layer was on the specimen surface. In the 3D analy- 
sis, the crack length was held constant (2c/w = 0.5) and no crack growth was considered. 
The minimum element size d at the crack front was 0.004 inches (d/c = 0.0016). As will be 
described later, the crack-front shapes for the non-straight cracks were determined from exper- 
imental observations. To model the nonstraight crack fronts, the mesh shown in Figure 4(b) 
was modified. The z-coordinates of a portion of the mesh surrounding the crack tip (0.7c < x 
< 1.3c and 0.0 < y < 0.3c) were changed to correspond to the measured crack fronts. 

RESULTS AND DISCUSSIONS 

Blunt-Notch Tests and Analyses 

Because the FEA codes employ the small-strain assumption, their ability to model large-scale 
plastic deformations was a concern. The blunt-notch specimen, Figure 2(b), was used to ex- 
amine this issue. In this specimen, the notch length remains fixed while the specimen under- 
goes intense plastic deformations. Displacements were measured at the center of the notch and 
at both notch roots during monotonic loading to failure. Figure 5 shows a comparison of ap- 
plied stress against notch-root displacement from the test (symbols) and from 2D FEA (solid 
curve). The -displacements were measured and calculated at the locations shown in the in- 
sert. The solid symbol is the last reading before the specimen failed. The analysis agreed 
quite well with the test results especially after net-section yielding (the plateau). In the early 
stages of deformation, the analysis was about 10 percent lower than the measured displace- 
ments indicating that the finite-element model was probably too stiff. Similar agreement be- 
tween measured and predicted centerline displacements was found for both the 2D analysis 
with a sharp notch and the 3D analysis with a non-growing crack. These results confirm that 
both the 2D and 3D analyses are able to model the large-scale plastic deformations that are 
present during the fracture of thin-sheet aluminum alloys. 

Stable Tearing Tests and Analyses 

The results of three M(T) fracture tests using three material thicknesses are shown in Figure 6 
where the applied stress (S) is plotted against crack extension. All of these specimens were 
fatigue precracked at a low stress level, then loaded monotonically until failure. These test re- 
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suits were used to determine the critical CTOA (vj/ c ). A value of \| / c was selected so that the 
FEA would produce a maximum load equal to the average of the maximum loads of these 
tests. The critical angle \|/ c was 6. 1° under plane-stress conditions. Calculations from the 
FEA using vy are shown by the solid curve in Figure 6. The agreement between the calcu- 

V 

lations and the test data near and beyond maximum load is reasonable considering that the net- 
section stresses S n are between the yield stress and the ultimate tensile strength. Calculations 
under plane-strain conditions using the same critical angle are shown by the dashed curve. 
Surprisingly, these results agreed well over the complete range of data. However, this agree- 
ment is fortuitous because the thin-sheet mat erial is definitely not in a condition of plane 
strain. 

Tests were conducted on specimens that were fatigue precracked at both high and low stress 
levels. The analyses and test data are compared in Figure 7 for the 0.07-inch thick material. 
The solid symbols show the measured crack extension results after the high precracking history 
(S max = 22.5 ksi). These results show that the higher fatigue stress delayed the crack initia- 
tion point to an applied stress level slightly higher than the maximum fatigue stress. The 2D 
FEA results for the high precrack stress, shown by the solid curve, agreed very well with the 
test data. In the analysis, the high fatigue stress caused fatigue crack closure to develop while 
the crack grew under fatigue. The crack tip was shielded by a wake of material that had been 
plastically deformed and, thus, a higher applied stress level was required to reach the critical 
angle vy c . The open symbols and the dashed line show the results for the low precrack stress. 
The test results and the analysis also agreed well for low fatigue stress levels. As previously 
discussed, a crack grown under a low fatigue stress will have to go through a transition to the 
shear mode fracture in the early stages of crack growth. Additionally, as will be discussed 
later, in the early stages of growth, the crack growth lags on the surface and grows more in the 
interior. Thus, the measured crack extensions on the surface may not be indicative of the av- 
erage crack extension. The high fatigue stress test did not have to transition to shear mode 
fracture and, as will be shown later, crack growth through the thickness is nearly uniform 
during shear mode fracture. 

Figure 8 shows a comparison of measured and predicted crack extensions for a smaller width 
M(T) specimen using vy c = 6. 1°. These results are quite similar to those shown for the wide 
specimens, that is the analysis (curve) overpredicted crack extension before maximum load. 
Again, crack tunneling may be the reason for the disagreement. However, the analysis pre- 
dicted a maximum applied stress within 4 percent of the experimental results. These results 
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indicate that the critical CTOA is able to predict the effects of specimen size on maximum 
failure loads. 

A comparison of measured and predicted load-line displacements (displacements remote from 
the crack) for the wide M(T) specimens is shown in Figure 9. The symbols show the meas- 
urements made on the three specimen thicknesses. Here the thinnest specimen (B = 0.05 
inches) tended to show higher displacements beyond maximum applied stress but this may 
have been due to the thinner specimen reaching a higher maximum applied stress level than the 
other specimens. Again, two analyses were made using the critical CTOA \|/ c of 6.1°. One 
analysis was made under plane-stress conditions (solid curve); the other was made under plane- 
strain conditions (dashed curve). As expected, the plane-strain analysis was too stiff. How- 
ever, beyond the maximum applied load, the plane-stress analysis predicted more displacement 
than was observed in the tests. The reason for this discrepancy is still under investigation. It 
may be related to constraint effects because the test data fell between the plane-stress and 
plane-strain results but closer to the plane-stress results. 

Comparison of Measured and Calculated CTOA 

As shown previously in Figure 2, the CTOA was measured during the fracture process on 
several M(T) specimens and one three-hole-crack specimen. The results for one of the M(T) 
specimens is shown in Figure 10. The critical angle v|/ c is plotted against crack extension, 
Ac. For each increment of crack extension, one to three values of \j/ c were measured at 
varying distances behind the crack tip. The vertical lines indicate the location of the initial 
crack tip Cj, the transition to shear mode fracture and the peak load. The measured V c be ‘ 
fore the transition are between 7.5 and 15° but after transition they appear to have an average 
value of about 6.5°. The solid horizontal line is the v|/ c value (6.1°) determined from the 2D 
FEA. The dashed horizontal line is the measured angle from the FRASTA method (20] on a 
thinner 2024-T3 aluminum alloy sheet. In the transition region, the agreement between the 
present measurements and the FRASTA method is quite good. The overall results from the 
present measurements are not vastly different from those presented in the literature for other 
materials. 

The measured vy values from a three-hole-crack specimen are shown in Figure 11. Again, 

v 

xj/ is plotted against crack extension. The vertical lines indicate the locations of the initial 

V 

crack tip Cj, the transition to shear mode fracture, the region of decreasing crack driving 
force (which begins at the location of the peak stress-intensity factor) and the center-line of the 
large holes. Each data point shown is the average of three to four \p c measurements made at 
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each value of crack extension. The trends are, again, similar to those shown in Figure 10. 
The critical vj/ determined from the analysis of the M(T) specimens is shown by the solid 

V” 

horizontal line at 6. 1°. Again, higher values of \j/ c were measured at crack initiation. The 
experimental and analytical results are in excellent agreement after the transition region. 

Crack Front Shape 

To understand the high CTO A values measured during crack initiation, a series of test were 
conducted to measure the crack front shapes during the fracture process, especially in the early 
stages before the transition to the shear mode fracture. Figure 12 shows measurements for 
crack-front shapes from three tests for specimens with a thickness B of 0.09 inches. The sym- 
bols represent the average values from the left and right crack tip measurements. The curves 
shown in the figure were used to model the crack fronts in the 3D FEA. The solid circles 
show the tunneling that would typically result from the fatigue precracking of the specimen. 
The specimen was then monotonically loaded and the specimen surface in the crack-tip region 
was observed using a high-power microscope. When a slight amount of crack extension was 
observed, the loading was stopped and a fatigue marker loading was applied. The specimen 
was then fractured to reveal the crack shape as shown by the open circles (Stable 1). A larger 
amount of crack extension was measured in the interior compared to the values at the free 
surfaces. A second specimen was tested where a greater amount of stable crack extension was 
allowed before fracturing the specimen; Stable 2 shows the crack shape that was measured. 
Even with the extreme tunneling found in this case, the crack front still lies within the flat 
fracture region. The cross-hatched region shows the development of the transition from the 
flat fracture to the shear mode fracture region. In the third test, the Stable 3 crack front was 
measured during the transition to the shear mode fracture; the crack front shape is becoming 
nearly straight through the thickness, although now on a 45° plane. These results suggest that 
in the early stages of stable crack growth, even in thin-sheet aluminum alloy, the fracture 
process is truly three dimensional. 

Three-Dimensional Finite-Element Analyses 

The 3D FEA code ZIP3D was used to analyze the influence of crack shape on crack-tip-open- 
ing displacement. Comparisons were made for four measured crack front shapes (Straight, 
Fatigue, Stable 1 and Stable 2) in the flat crack region, as shown in Figure 12. For the 3D 
analyses, the crack length was held constant for all crack front shapes and no crack extension 
was modeled. Crack-tip-opening displacements were compared at a distance d = 0.004 inches 
behind the crack tip; this distance corresponds to one element length (d/c = 0.0016). Normal- 
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ized elastic displacements for the four crack shapes are shown in Figure 13, where E is 
Young's modulus, S is applied stress, B is the specimen thickness (B = 0.09 inches), and V is 
the crack-tip-opening displacement. These results show that even modeling the small amount 
of tunneling found in the fatigue crack front will greatly change the crack mouth profile. All 
three of the curved crack fronts show a significantly larger displacement on the specimen sur- 
face (z/B = 0.5) that in the center of the specimen thickness (z/B = 0.0). 

Experimentally, stable crack growth was observed at the surface at an applied stress of about 
15 ksi when the crack front shape was somewhere between the Fatigue crack shape and the 
Stable 1 crack shape. The 3D elastic-plastic analysis was used to calculate the crack-tip-open- 
ing displacements corresponding to both the Fatigue and Stable 1 crack fronts at an applied 
stress of 15 ksi, as shown in Figure 14. The two curves are similar. The displacements calcu- 
lated for the Fatigue crack front are 31 % greater on the surface compared to the values at the 
center; for the Stable 1 crack front, the displacements are 44% greater on the surface com- 
pared to the center. If the 6. 1° determined earlier is assumed to govern the crack growth in 
the interior, then the surface CTOA should be 8 to 9°. This agrees with the trends in vj/ c 
measured prior to the transition, shown earlier in Figures 10 and 11. Higher values of v|/ c 
were measured on the surface but lower values in the interior would control the crack growth. 
Crack-tip-opening displacements were also calculated for the Stable 2 crack shape at a higher 
applied stress level (S = 27 ksi), as shown in Figure 14. Although these crack-tip displace- 
ments show variation through the thickness, it appears that with the changing crack front 
shape, the displacements are becoming more uniform through the thickness. 

CONCLUDING REMARKS 

A two-dimensional, elastic-plastic, finite-element analysis was used with a critical crack-tip- 
opening angle (CTOA) fracture criterion to model stable crack growth in thin-sheet 2024-T3 
aluminum alloy under monotonic loading after precracking at either low or high cyclic stress 
levels. Tests were conducted on three types of specimens: middle-crack, three-hole-crack and 
blunt-notch specimens. All specimens were subjected to remote tensile loading and tested un- 
der displacement-control conditions. The specimen thicknesses ranged from 0.05 to 0.09 
inches. An experimental technique was developed to measure CTOA during crack-growth in- 
itiation and stable tearing using a high-resolution video camera and recorder. Crack-front 
shapes were also measured during initiation and stable tearing by using a fatigue marker-load 
technique. Three-dimensional elastic-plastic finite-element analyses of these crack shapes for 
stationary cracks were also conducted to study the crack-front opening displacements. 
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Predicted load against load-line and notch-tip displacements under plane-stress conditions for 
the blunt-notch specimen agreed well with test results for large-scale plastic deformations. 
Comparisons made between measured and calculated load against crack extension behavior on 
middle-crack tension specimens agreed well even as the net-section stresses approached the ul- 
timate tensile strength. The analyses were able to predict the effects of specimen size and 
precracking stress history on stable tearing. Predicted load against load-line displacements for 
middle-crack specimens agreed well with test results up to maximum load but the analyses 
tended to overpredict displacements beyond maximum load in displacement-controlled tests. 
During the initiation phase, the measured CTOA values from surface observations were high 
but they became nearly constant after a small amount of stable tearing. The average measured 
value of CTOA agreed well with the calculated value from the finite-element analysis. The 
high CTOA values measured at the sheet surface during the initiation phase may be associated 
with the crack tunneling that was observ ed i n the tests. Three-dimensional analyses for non- 
straight crack fronts gave much higher displacements near the free surface than in the interior. 
The analysis methodology presented here may be used in global-local structural analysis 
methods to predict stable crack growth in cracked structural components. 
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(a) Middle-crack (b) Blunt-notch (c) Three-hole-crack 

Figure 1. - Specimen configurations tested and analyzed. 
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Figure 2. - Computer-enhanced high-resolution video frame of stably tearing crack. 


14 





Fatigue surface 



Figure 3. - Schematic of fracture surface indicating transition from a flat 
to a slant crack plane. 



Elements = 4185 
Nodes = 2287 
d = 0.01875 In. 
w = 5.9 In. 


********** 

********** 

********** 

********** 

********** 



Layers = 6 
Elements = 5706 
Nodes = 7203 
d = 0.004 In. 
w = 5 In. 



(a) Two-dimensional model (b) Three-dimensional model 

Figure 4. - Typical finite-element models of middle-crack specimens. 
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Figure 5. - Comparison of measured and predicted notch-tip displacements for 
blunt-notch specimen statically pulled to failure. 



Figure 6. - Comparison of measured and predicted stable crack growth in wide 
M(T) specimens using critical CTOA criterion. 
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Figure 7. - Comparison of measured and predicted stable crack growth in wide 
M(T) specimens for low and high precracking stress. 



Figure 8. - Comparison of measured and predicted stable crack growth in a 
small width M(T) specimen using critical CTOA criterion. 
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Figure 9. - Comparison of measured and predicted load-line displacements in 
wide M(T) specimens using critical CTOA criterion. 


r FEA (Present study) 

-- Ref. 20 (B = 0.04 in.) 


15 


Tp 

c 10 h 

deg. 


d*o § 


Peak load 


8 


2024-T3 

Middle-crack tension 
w = 5.9 in. 

Cj = 2 in. 

B = 0.09 in. 


% o 


~u 


O O ^ 
CL ° 


I H 

C: Transition region 


o 


t r 


0 

-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Crack extension, Ac, in. 


Figure 10. - Comparison of measured and calculated CTOA from M(T) specimens. 
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Figure 11.- Comparison of measured CTOA from a three-hole-crack specimen 
with calculated CTOA from M(T) specimens. 


, Straight Crack Front 


Shear Fracture 
(45°) 



Flat 

Fracture 


Drawn to Scale 


Figure 12. - Crack-front shapes measured during fatigue crack growth and 
stable tearing for M(T) specimens. 
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Figure 13. - Effect of crack front shape on calculated elastic crack-opening 
displacements for M(T) specimens using 3D FEA. 



Figure 14. - Effect of applied stress on elastic-plastic crack-opening 
displacements for M(T) specimens using 3D FEA. 
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